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We study the potential between static SU(4) sources using the Model of Thick Center Vortices. 
Such vortices are characterized by the center elements zi — i and Z2 — zf. Fitting the ratios of 
string tensions to those obtained in Monte-Carlo calculations of lattice QCD we get /2 > fi, where 
/„ is the probability that a vortex of type n is piercing a plaquette. Because of 22 = zf vortices of 
type two are overlapping vortices of type one. Therefore, /2 > /i corresponds to the existence of 
an attractive force between vortices of type one. 



I. INTRODUCTION 



It is a well known fact that Quantum Cliromodynamics (QCD) is truly described by the SU(3) gauge group. How- 
ever, the mechanism of confinement is still under intensive discussion and is an attractive open subject. Studying the 
SU(N) gauge theories with N greater than three, helps us to better understand different aspects of QCD. Confinement 
implies that the potential between static sources has a linear term 

V~crR (1.1) 

where a is the string tension and R is the distance between the two sources. SU(N) gauge theories with > 3 
have some characteristics that SU(3) does not have. For example, there exists only one universal string tension for 
the SU(3) gauge theory but for SU(N) with iV > 3, the number of independent stable string tensions are equal to 
int[^]. The values of these string tensions will constrain the details of the confinement models. Both meta-stable 
string tensions which are the string tensions at intermediate distances and the stable string tensions which are the 
string tensions at large distances, are important, in this respect. Stable string tensions depend on the N-ality k 
of the corresponding representation. In fact, at large distances where the distance between two sources increases, 
a pair of gluons pops out of the vacuum and couples with the initial sources. They do not change the N-ality of 
the representation but reduce the dimension of the representation to the lowest dimension of the corresponding N- 
ality. Therefore, stable (asymptotic) string tensions depend on the N-ality of the representation. On the other hand, 
at intermediate distances, the string tensions depend on the dimension of the representation. For this region, the 
potential energy is not large enough to produce a pair of gluons and therefore, the dimension of the original sources 
does not change. 

For stable strings, there are some existing theories describing the ratio of ^ where erf and ak are the string tensions 
for fundamental quarks and for quarks of representations with N-ality k, respectively. The linear potential between 
two static sources may be explained by forming a chromoelectric flux tube carrying charge in the center Zj^ of the 
gauge group. The most trivial idea is that the total flux is carried by k independent fundamental tubes. Then 

<7k = ka{ ; fc min{fc, N - k} (1.2) 

Because of charge conjugation we get ak = ON-k- Thus, for the SU(3) gauge group, 02 = <Ji and one universal string 
tension is obtained. If one wants to study theories with more than one string tension, one has to study SU(N) with 
> 3. The asymptotic Casimir scaling is another theory which claims that [i| 

au ^ k{N~k) 

CTf iV - 1 ^ ' ^ 

Calculations in brane M-theory predict Sine-law scaling 



(1.4) 



On the other hand, for intermediate distances, lattice calculations show that the string tensions are roughly pro- 
portional to the eigenvalues of the quadratic Casimir operators 0] . In ref . ^ this phenomenon was dubbed "Casimir 
scaling" . The Casimir scaling regime is expected to extend roughly from the onset of confinement to the onset of 
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screening |5| . There is another argument about the linear part of the potential at intermediate distances which claims 
that the string tension in this region is proportional to the number of fundamental flux tubes embedded into the 
representation . The fundamental flux or string is the one that connects a fundamental heavy quark with an anti- 
quark. This idea is called "flux tube counting" . In general, the fundamental strings do not interact at very large N. 
Thus if one interprets the meta stable string tensions between the higher representation sources to be proportional to 
the number of flux tubes with jj corrections; then for large enough distances these meta-stable strings will decay into 
stable strings with given A^-ality, whose tension is likely to be described by the Sine formula or by Casimir scaling in 
some approximation which is not yet known so far. 

In this article we study the string tensions of SU(4) static sources, within the model of thick center vortices. In our 
previous calculations 0, we have shown that using only the first non-trivial type of the vortices of the SU(4) gauge 
group, two different asymptotic string tensions may be obtained, one for fc = 1 and another one for fc = 2. Their 
ratio approaches asymptotically ~ 2 in agreement with flux tube counting, see Eq. (|1.2p . If one wants to get 
other asymptotic string ratios, e.g. those corresponding to lattice calculations, one has to use both vortex types of 
the SU(4) gauge group. In the following sections, after giving a brief review of the role of vortices for the confinement 
of quarks, potentials between static sources are calculated using both types of SU(4) vortices. Because of the relation 
between the ratio of the stable string tension, where crt is the fundamental string tension, and the probabilities 

/i and /2 of piercing plaquettes by vortices of type one and two, one can determine the ratio of the probabilities, 
by fixing — from lattice calculations or from the theories related to Eqs. (|1.2p - (|1.4[) . Then, the induced potentials 
may be determined from the model. We show that the meta-stable string tension ratio ^ is almost independent 
of the ratio of the asymptotic string tensions, where ar is the string tension of representation r at intermediate 
distances. Approximate agreement with both Casimir scaling and flux tube counting is observed for all potentials at 
intermediate distances. The effect of the second vortex type which modifies the concavity of the potentials and the 
type of the interaction between vortices are also discussed. 



II. POTENTIALS FOR TWO TYPES OF THICK CENTER VORTICES 



The vortex model of QCD assumes that the vacuum of quantum chromodynamic is filled with vortices of finite 
thickness which carry magnetic fluxes corresponding to the center of the gauge group. In order that vortices have 
finite energy per unit length, their gauge potential at large transverse distances must be a pure gauge. The non-trivial 
nature of gauge transformations for producing the gauge potentials causes that the vortex cores have non-zero energy 
and makes vortices topologically stable. The potential energy between static sources induced by the vortices is Q 



N-l 



V{R) 



(2.1) 



xeA 



X is the location of the center of the vortex and C indicates the Wilson loop. A is the minimal area of the Wilson 
loop and the sum over positions x runs over all plaquettes in the plane of the loop. /„ represents the probability that 
any given unit area is pierced by a vortex of type n. which gives the information about the flux distribution and 
the contribution that a vortex with its center in a specific plaquette may have to the Wilson loop is given by 



Qr[a\ 



1 

dr 



Trexp[ia.i/''], 



(2.2) 



where dr is the dimension of the representation r and {i/,f , i = 1,2, — 1} are the generators in this representation 
spanning the Cartan subalgebra. The vector adx) depends on the fraction of the vortex core which is enclosed by 
the loop, thus on the color structure and the position of the vortex and the shape of the loop. For each SU(N) gauge 
group, there are A^ — 1 types of vortices, corresponding to the non-trivial center elements 2:„. Vortices of type n and 
N — n differ in the direction of the magnetic flux only and have the same probability /„. Thus 



fn = In-t 



and 



(2.3) 



Hence, among the three vortices of the SU(4) gauge group, two of them may be considered to be the same. Therefore, 

/i is equal to /a and Re5r[Q?|^'' (a;)] = KeGriaQ^x)], where the upper index (n) in a[?''(a;) indicates the type of the 
vortex. 

According to Eq. (|2.ip the induced potential between SU(4) sources may be written as 



V{R) = - ^ In {l - 2/i (l - ReGr [4'^ (2^)] ) - /2 (l - Re^, [a^ (x)] ) } 



(2.4) 



3 



Since /i = /a, the first term in equation (|2.4p is multiplied by 2. To determine Gr [a] from equation (|2.2p . iJ's and a's 
for the SU(4) gauge group should be determined. The defining (fundamental) representation of the three generators 
Hi of the Cartan sub-algebra may be chosen as 

i^i- ^(1,-1,0,0), 

^2 = ^(1,1,-2,0), (2.5) 
i/3 = ^(1,1, 1,-3). 

A center vortex completely linked to a Wilson loop, in the fundamental representation of the SU(N) gauge group, 
has the effect of multiplying the Wilson loop by a non-trivial center element z„ — exp{^^} 

W{C)^eM^}W{C) n-1,2,... ,7V-1. (2.6) 

For the group SU(4), a vortex of type n = 1 requires zi = i and a vortex of type n = 2 requests Z2 = —1. The 
percentage of linking we describe by a function r]{x) which is zero when the vortex does not touch the loop and equal 
to one if the vortex is entirely contained within the loop. Then, the influence of a vortex of type n on the Wilson loop 
is 

^^[a(»)] = lTrexp{ia(")(.x)i?}, with a^"\x) = a^^ler]{x) (2.7) 

A possible choice for the unit vector e and amix is 

e= (0,0,1), 41), = ^V6, a(2)^ = 27rV6, (2.8) 

leading to 

= ^Tr exp{iri{x)mTVQH3} = ^ [3 exp{i7y(x)n^} -f exp{i?7(a;)n^}] (2.9) 

We can choose other unit vectors which lead to a permutation of the matrix elements of H3 or to a sign change of H3 
and leave the sum in Eq. (|2.9p unchanged, e.g. 

1^,1) with ei? = ^diag(3, -1,-1,-1) (2.10) 
We have used the following profile, see ref. 

rj{x) ^ ^[1 - ta.nh{ay{x) + ^)], (2.11) 

is the distance of the vortex center to the nearest timelike side of the loop. The parameter a is of the order of the 
inverse of the vortex core thickness [3] and the parameter b takes into account that vortices can not be completely 
linked to small Wilson loops. 

Therefore, the vectors o?^") (x) satisfy the following conditions 

1. Vortices which pierce the plane far outside the loop do not affect the loop. It means for fixed R, as x ^ 00, 

2. If the vortex core is entirely contained within the Wilson loop, then for vortex type one, |cf'^-'| — VGtt and for 
vortex type two, |q?(^)| = 2^/6it. 

3. As i?^ then 0. 

So far, we have determined the group valued functions Gr [o?] of equation (|2.2p . However, the two parameters /i and 
/2 are not specified yet. In the next section, we show that the asymptotic string tensions depend on the 4-ality class 
k of the Wilson loop and deduce from the asymptotic string tension ratios the probability ratios 4^. 



where 
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III. ASYMPTOTIC STRING TENSIONS 



The asymptotic string tensions can be determined from very large Wilson loops. In this regime, we can neglect 
the finite thickness of vortices and assume that vortices piercing the minimal area of the loop would insert a center 
element somewhere in the product of link variables. For the SU(4) gauge group, there are four center elements 

zo = 1 zi = exp (y ) Z2 = exp (ttI) 23 = 2*= exp (^)- (3.1) 

Which center element contributes depends on the type n of the vortex and the N-ality k of the representation r of 
the loop. The Young-tableaux of the lowest representations of 4-ality k for SU(4) and their dimensions are depicted 
in Fig. [TJ One can easily verify in this figure that the N-ality k modulo N is given by the number of quarks which is 

k = l fc = 2 fc = -l k = 
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FIG. 1: The Young tableaus for SU(4) representations of various dimensions and their N-ality k is shown up to four-quark 
states. 

equal to the number of squares in the Young tableau which is necessary to get such a representation. A representation 

k(r) . 

r contributes with a factor Qr — Zn ■ Hence, from equation (|2.ip . the asymptotic potential for a representation r of 
N-ality k reads 



^ In |l - ^ - Re{z'^^(^y)]\ = eM-Aar}. 

xeA I n=l J 



Vr{R) 

Using /i — from Eq. (|2.3p we get for the fundamental representation of SU(4) 



<Jf = - ln{l - 2/1 [1 - Rc(zi)] - /2 [1 - Rc(z2)]} = 

- - In [1 - 2/1 - 2/2] ~ 2/1 + 2/2 for (A , /2 « 1) 

For the diquark sources with 4-ality k — 2, like the sextet and the decuplet follows 

CT6 = fiio = - In {1 - 2/1 [1 - Re{zf)] - ^ [l - Re(z2')] } = 
= - In (1 - 4/1) ^ 4/1 for (/i , /2 « 1) 



(3.2) 



(3.3) 



(3.4) 
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Finally, we see that the asymptotic string tension of the 4-ality fc = representations, which can be build by sources 
with the same numbers of quarks and antiquaries, is equal to zero. The most important member of this class is the 
adjoint representation 

aadj = - In { 1 - 2/1 [1 - Re(zO)] - /a [l - Re(zO)] } = 0. (3.5) 
To summarize, the SU(4) asymptotic string tensions behave for small piercing probabilities like 

(7f = 2/i + 2/2, a6 = 4/i, aadj=0. (3.6) 

This leads to the ratio 



C^f /l + /2 



(3.7) 



of the asymptotic string tension for representations with 4-ality fc = 2 and 4-ality fc = 1. 

This behavior that the asymptotic string tensions are equal for all quark sources of the same N-ality one can picture 
as pairs of gluons which pop out of the vacuum and screen the original sources. Therefore, the SU(4) gauge group 
has two universal asymptotic string tensions and Eq. p.7p gives the ratio of these string tensions within the model 
of thick center vortices in terms of /i and /2, of the probabilities of piercing plaqucttes by vortices of type n = \ and 
n = 2. 

In the next section, we use the ratio of ^ from lattice calculations and from the predictions (|1.2p - (|1.4p to fix the 
ratio of by equation p.7p . Then the potentials between static sources are calculated and discussed. 



IV. RESULTS AND DISCUSSION 



In our previous calculations fT^I for the gauge group SU(4), we assumed that vortices of type n = \ pierce plaquettes 
with the probability /i and no vortices of type n = 2 are present, /2 = 0. From Eq. (j3.7p we get under these 
assumptions erg — laf . This agrees with Eq. ()1.2|) and corresponds to the picture that a string of N-ality A; = 2 is built 
from two non-interacting fundamental strings. Lattice calculations do not confirm this scenario, since they predict 
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10|: —- 1.370(20) — > ^= 0.0460(7), 

L. Del Debbio et. al. [9]: — ~ 1.403(15) — > — = 0.0426(5). 

CTf /i 

On the other hand, the limit /i — f2 leads to erg — ; in obvious contradiction to the lattice results. 

We fix now the ratio of to the values (|4.ip predicted from the lattice data and adjust the absolute values of the 
probabilities /i, /2 ai, 02, 61 and 62 such that the potentials at intermediate distances become linear and get the 
best agreement with the intermediate string tensions of the lattice data. For both vortex types, the general form of 
the vortex profile introduced in equation (|2.1ip has been used. Fig. [5] shows the potentials of various representations 
versus R in the range of i? e [1, 100], using 

/i = 0.1, /2 = 0.046, ai = 0.05, h = 4.0, 02 = 0.025, 62 = 8.0. (4.2) 

Indices one and two refer to vortices of type one and two. As one can see in Fig. [5] at large distances, zero 4-ality 
representations, like 15 (adjoint), 20n and 35 are screened. For non-zero 4-ality representations we get two different 
asymptotic string tensions. The potentials between the sources with dimension 20^ and 20nis become parallel to that 
of the fundamental representation and the potentials of the diquark representations (6 and 10) get the same slope. 

As expected, we find for all the potentials in Fig. [2] a linear region at intermediate distances. Table [J shows the 
ratios — for the corresponding intermediate string tensions cr^ of some of these representations for /i ~ 0.1 and for 
some values of /2 . The errors in parentheses ^ show the systematic errors due to slight changes of the linear regime. 
The errors in square brackets are due to the uncertainties of the lattice data we have used. The first row contains 
the results of our previous work [7] with /2 = 0, the second and the third line the above discussed results from the 
comparison with the lattice calculations of refs. Q and and the forth row the case /i = /2. It is interesting to 
compare these values with the absence of vortices of type n = 1 in the fifth row (/i = and /2 = 0.1) and with the 
fits for the asymptotic string tension to the predictions of the Casimir scaling law in Eq. (|1.3p and the Sine law scaling 
in Eq. (jl.4p in lines six and seven. It is clearly seen that slight changes of the piercing probabilities /2 have only a 
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ratios 


£6 


<^ 15 
<Tf 


^ 10 


°'20a 




^ (asymptotic) 


fl = 0.1, /2 = 


1.50(5) 


1.73(6) 


1.88(5) 


2.65(10) 


3.21(27) 


2 


fl = 0.1, /2 = 0.0426[5] 


1.51(11)[2] 


1.71(9)[2] 


1.88(10)[2] 


2.57(15)[3] 


3.22(18)[3] 


1.403[15] 


fl = 0.1, /2 = 0.0460[7] 


1.50(8)[2] 


1.69(14)[3] 


1.86(8) [3] 


2.52(10)[4] 


3.16(14)[5] 


1.370 [20] 


/l = 0.1, /2 = 0.1 


1.53(15) 


1.71(16) 


2.00(15) 


2.72(21) 


3.39(26) 


1 


/l =0, /2 = 0.1 


1.46(19) 


1.48(9) 


1.67(12) 


2.18(17) 


2.64(18) 





/i = 0.1, /2 = 0.050 


1.43(20) 


1.71(12) 


1.88(13) 


2.61(18) 


3.22(23) 


1.333 


/i = 0.1, /2 = 0.042 


1.37(20) 


1.78(18) 


1.99(19) 


2.71(26) 


3.34(33) 


1.414 


Casimir ratio 


1.33 


2.13 


2.4 


4.2 


6.4 




no. of fund, fluxes 


2 


2 


2 


3 


4 





TABLE I: The string tension ratios ^ for some SU(4) representations r at intermediate distances are shown in the flrst five 
columns and the values of ^ in the asymptotic region in the last column. The first four rows treat the potentials for the 
piercing probabihties fi =0.1 and increasing values of /2, for our previous work [7| in the first row, for the fits to the lattice 
calculations of refs. and[l3l in the second and third row and for the case /i = /2 in the forth row. Line five shows the ratios 
in the absence of vortices of type n = 1, line six fits to the predictions of the Casimir scaling law in Eq. (|1.3|l and line seven fits 
to the Sine law scaling in Eq. (|1.4p in the asymptotic region. These values can be compared with the the ratio of eigenvalues ^ 
of the quadratic Casimir operator i n line eight and with the number of fundamental fiuxes in line nine. The errors in square 
brackets are due to the uncertainties of the lattice data we have used, the errors in parentheses ^ show the systematic errors 
due to slight changes of the linear regime. 

weak influence on the intermediate string tensions. This situation at intermediate distances differs drastically from 
the behavior at asymptotic distances, which was discussed above and is shown for comparison in the last column of 
Tab. [H The asymptotic string tensions of the lattice calculations in refs. Q and [l3| can be reproduced only with 
an appropriate density of vortices of type n — 2. The ratio of eigenvalues ^ of the quadratic Casimir operator and 
finall y th e number of fundamental fluxes are represented in lines eight and nine. The lattice results in Table 17 of 
ref. [lO| deviate for the representations 6 and 10 (there is no report on higher dimensional representations) by about 
3.5% and 15% from Casimir scaling, our results in lines 2 and 3 of Table [J by about 14% and 23%. The agreement 
with Casimir scaling gets worse as the dimensions of the representations increase. This fact is also observed in ref. [l^ 
for the representations 6 and 10. Considering that the thick-center-vortices model does not reproduce the Coulombic 
part, the behavior of the linear part we get from the model is satisfying. 

For values of /i and /2 which are small enough compared to one, we can expand the logarithm in Eq. (j2.4p around 
one and separate the contributions for the potentials produced by the two types of vortices. For our choice of the 
piercing probabilities in Fig. El /i = 0.1 and /2 = 0.046, the violation of this additivity of the two contributions to 
the potential by the exact expression (|2.4p is invisible, see Fig. [31 which compares the full decuplet potential with 
the potential produced by vortices of type n = 1 (/i = 0.1 and /2 = 0.0) and vortices of type n — 2 {fi — 0.0 and 
/2 = 0.046). It is nice to see that the vortices of type n — 2 produce k — 2 potentials which are asymptotically 
screened as predicted by Eq. (13. 4|) . This can be understood by the QCD analogon of the Aharanov-Bohm effect. 
Carrying a two-quark source around a vortex corresponding to a Z2 = = — 1 color magnetic flux leads to a phase 
of z| = 1 , to a screened potential. It is easily understandable from Eq. (|2.4p that the value of this screened potential 
is roughly proportional to the piercing probability /2. Due to this screening effect, vortices of type n — 2 increase the 
string tension at intermediate distances only. 

There is a probability that non-interacting vortices of type n — 1 have the same position and due to Eq. (|2.8p 
are identified as vortices of type n = 2, /2 = /f. A comparison to the lattice calculations, see Eq. (|4.1[) . leads to 
f2 > fl - This indicates an attraction between parallel vortex fluxes. 

V. CONCLUSION 

Confinement is one of the most interesting features of QCD which has been studied by both lattice gauge theory 
and phenomenological models. The model of thick center vortices is one of the phenomenological models which has 
been fairly successful to explain the linear part of the potentials. For the SU(N) gauge groups with > 4 there exist 
vortices with different quantized fluxes. In this article we have studied the gauge group SU(4) which has two types 
of vortices, vortices of type n = 1 with magnetic flux corresponding to the first non-trivial center element zi — i 
of SU(4) and vortices of type n = 2 with a flux corresponding to Z2 — —1. We have shown that the ratio /2//1 of 
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FIG. 2: Potentials for various representations of SU(4) using both types of vortices for piercing probabilities /i = 0.1 and 
/2 = 0.046, see Eq. (|4.ip . The dimensions of the representations are indicated. 



the probabilities for the piercing of plaquettes by vortices of both types determines the ratios of asymptotic string 
tensions of 4-ahties A; = 0, 1 and 2. We underhne that the lattice results for the ratios of these string tension can be 
explained only by using both types of vortices . Adjusting these ratios to the slightly different results of refs. 0] and 
[10| . the general features of the potentials and the string tensions at intermediate distances are yet indistinguishable. 

Using vortices of type n — 2 only, the six and ten dimensional representations, which have 4-ality k — 1 and are 
possible two-quarks states, are screened at large distances. This results from the multiplication of the Wilson loop 
holonomy in a two-quark state by z| = 1. More generally, it follows that the k = 1 asymptotic string tension is 
independent of the probability /2. The value of fi influences the fc = ±1 asymptotic string tensions only. 

Because of Z2 — z\q. vortex of type n = 2 corresponds to two overlapping vortices of type n = 1. The analysis of the 
probabilities /i and gives an information about the interaction of vortices. Non-interacting vortices of type n = 1 
are described by fi = , whereas /2 > indicates vortex attraction and /2 < vortex repulsion. A comparison of 
our results with Monte-Carlo calculations indicates that vortices attract each other. 

A consideration of vortices of type n = 2 modifies the concavity of the potentials which has been observed in 
previous calculations of SU(2), SU(3) and SU(4) potentials with only one type of vortices, see refs. 0,01,111 and 
which is not physical [Tlj . We conjecture that closed unquantised random magnetic flux lines of rather small size 
allow to remove this concavity and to introduce a Coulombic contribution with the correct sign. 
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FIG. 3; For piercing probabilities /i and f-z small compared to one the contributions of the two vortex types of SU(4) are 
almost completely additive as one can see in these decuplet potentials for the indicated piercing probabilities. Remarkable is the 
screening of potentials of 4-ality k — 2 produced by vortices of type n = 2 (/i = 0.0 and fz = 0.046) which is understandable 
by a QCD-Aharanov-Bohm effect. 
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